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Abstract 

We define truncated Mellin moments of parton distributions by restricting the integration range 
over the Bjorken variable to the experimentally accessible subset xq < x < 1 of the allowed 
kinematic range < a; < 1. We derive the evolution equations satisfied by truncated moments 
in the general (singlet) case in terms of an infinite triangular matrix of anomalous dimensions 
which couple each truncated moment to all higher moments with orders differing by integers. 
We show that the evolution of any moment can be determined to arbitrarily good accuracy by 
truncating the system of coupled moments to a sufficiently large but finite size, and show how 
the equations can be solved in a way suitable for numerical applications. We discuss in detail 
the accuracy of the method in view of applications to precision phenomenology. 
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1 The use of truncated moments 



The needs of accurate phenomenology at current and future hadron colhders have recently led 
to the development of more refined tools for the QCD analysis of collider processes [|l|. A 
detailed understanding of parton distributions and their scaling violations is in particular an 
essential ingredient of such phenomenology [|l], Scaling violations of parton distributions 
are described by renormalization group equations for matrix elements of leading twist opera- 
tors, whose anomalous dimensions are currently fully known to next-to-leading order [Q, and 
to next-to- next-to-leading order for a handful of operators @]. Moments of parton distribu- 
tions are related to moments of deep-inelastic structure functions by Wilson coefficients, which 
have also been computed up to next-to-next-to- leading order |^. Moments of structure func- 
tions, however, cannot be measured even indirectly, since they are defined as integrals over 
all kinematically allowed values of x, and thus require knowledge of the structure function for 
arbitrarily small x, i.e. arbitrarily large energy. 

There is of course a well-known solution to this problem, which consists of using the Altarelli- 
Parisi equation p to evolve parton distributions directly: the scale dependence of any parton 
distribution at xo is then determined by knowledge of parton distributions for all x > xq, i.e., 
parton evolution is causal. In fact, through a judicious choice of factorization scheme |0, || all 
parton distributions can be identified with physical observables, and it is then possible to use 
the Altarelli-Parisi equations to express the scaling violations of structure functions entirely in 
terms of physically observable quantities. It is, however, hard to measure local scaling violations 
of structure functions in all the relevant processes: in practice, a detailed comparison with the 
data requires the solution of the evolution equations 

What is usually done instead is to introduce a parametrization of parton distributions, 
and then solve the evolution equations in terms of this parametrization. The idea is that a 
parametrization fitted to the data will reproduce them in the experimentally accessible region, 
but it will also provide an extrapolation, so that the evolution equations can be solved easily, 
for instance taking Mellin moments. The results in the measured region should be independent 
of this extrapolation since, by the Altarelli-Parisi equation, measured scaling violations are 
independent of it. It has however become increasingly clear that in practice this procedure 
introduces a potentially large theoretical bias, whose exact size is very hard to assess [0, p!0| . 
First, the very fact of adopting a specific functional form constrains not only the extrapolation 
in the unmeasured region, but also the allowed behavior at the boundary of the measured 
region. Especially when data are not very precise, it can be seen explicitly [11| that rather 
different results are obtained simply by changing the functional form used to parametrize parton 
distributions. Furthermore, it is very hard to assess the uncertainty on the best-fit functional 
form of the parton distributions, essentially because of the very nonlinear and indirect relation 
between the data and the quantity which is parametrized. Hence, the need to go through such 
a parametrization makes it very hard to assess the uncertainty on the desired result. 

Various methods to overcome these problem have been discussed in the literature. One 
possibility is to minimze the bias introduced by the parton parametrization, by projecting 
parton distributions on an optimized basis of functions, such as suitable families of orthogonal 
polynomials P]. A more ambitious proposal is to construct the probability functional for parton 
distributions through Bayesian inference applied on a Monte Carlo sampling of the relevant 
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space of functions ||T0[- The probability functional then summarizes in an unbiased way all the 
available experimental information, and can be used to determine the mean value and error 
on any physical observable. For many applications, however, there is a simpler, although more 
limited option, which consists of dealing directly with the experimentally accessible quantities. 
Indeed, consider two typical problems in the study of scaling violations: the determination of 
and the determination of a moment of the gluon distribution, for instance the first moment 



of the polarized gluon distribution [|TT|, which gives the gluon spin fraction. Manifestly, in 
the latter case only the contribution to the moment from the measured region Xq < x < 1 
is accessible experimentally, and it would be useful to be able to separate in a clean way 
the measured quantity from the extrapolation to the unmeasured region, which is necessarily 
based on assumptions. It is then natural to study the scaling violation of these measurable 
contributions to Mellin moments, i.e. truncated moments. Similarly, as can be determined 
from the evolution equation of truncated moments |]T^, without any reference to the behavior 



of the structure function in the region in which it is not measured, and in principle without the 
need to resort to a specific functional parametrization. 

It turns out [|12| that scaling violations of truncated moments are described by a triangu- 
lar matrix of anomalous dimensions which couple the n-th truncated moment to all truncated 
moments of order n + k, where k runs over positive integers. This means that truncated mo- 
ments share with full moments the property that evolution equations are ordinary first-order 
differential equations, and not integro-differential equations, as for the parton distributions 
themselves. Unlike full moments, however, truncated moments can be measured without ex- 
trapolations. The price to pay for this is that their evolution equations do not decouple (unlike 
the case of full moments); however, the causal nature of the evolution implies that the evo- 
lution of each moment is only affected by higher order moments. Furthermore, the series of 
couplings to higher moments converges, and thus it can be truncated to any desired accuracy. 
The problem then reduces to the solution of a system of ordinary differential equations coupled 
by a triangular matrix, whose initial conditions are (measurable) truncated moments. 

This gives a simple solution to both problems mentioned above: and the truncated mo- 
ment of the gluon can be determined directly from the observed scaling violations, without 
having to go through an intermediate parton parametrization. A model independent error 
analysis can then be performed, provided only that data for the truncated moments and their 
errors are available. Of course, these could be extracted directly by summing over experimental 
bins, if a sufficiently abundant data set were available. In practice, however, it is more con- 
venient to manipulate a smooth interpolation of the data. It turns out to be possible to do 
this without invoking an explicit functional parametrization for the measured structure func- 
tions, by constructing neural networks which are trained to simulate all available experimental 
information, including statistical and systematic errors and correlations. By means of these 
neural networks it is then easy to compute the observed truncated moments and their errors, 
which can be further used to perform a phenomenological analisys of scaling violations, free of 
theoretical bias. The application of the method of neural networks to the parametrization of 
structure functions is currently under way and will be presented in a separate publication. 

Evolution equations for truncated moments were presented in ref. in the simplest (non- 
singlet) case, along with a preliminary study of the viability of the method. It is the purpose of 
this work to present a full treatment of the method of truncated moments, suitable for future 
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phenomeno logical applications. In particular, in sect. 2 we will derive evolution equations for 
truncated moments in the general (singlet) case, and we will discuss their solution in a form 
which is suitable for numerical implementation. In sect. 3 we will consider numerical solutions 
of the evolution equations with typical quark and gluon distributions, and discuss the accu- 
racy of the truncation of the infinite system of coupled evolution equations. Problems and 
techniques of interest for future phenomenological applications are discussed in sect. 4. The 
appendices collect various technical results which are needed for the actual implementation of 
the methods discussed here in an analysis code: in particular, we give explicit expressions for 
the NLO singlet splitting functions in the DIS scheme, and we list all the integrals needed to 
compute their truncated moments. 



2 Evolution equations for truncated moments and their 
solutions 

The dependence of parton distributions Q'^) is governed by the Altarelli-Parisi (AP) 
equations 

Q^) = f K asiQ')] Q') , (2.1) 

at ziT Jx y \y J 

where t = logQ^/A^. In the non-singlet case, q is simply one of the flavour non-singlet combi- 
nations of quark distributions, and P the corresponding splitting function. In the singlet case, 
g is a vector, whose components are the flavour-singlet combination of quark distributions, 

E{x,Q') = Y.q.{x,Q') (2.2) 

i=l 

and the gluon distribution g{x, Q"^). Correspondingly, P is a 2 x 2 matrix of splitting functions, 
given as an expansion in powers of ag- 

As is well known, upon taking ordinary Mellin moments convolutions turn into ordinary 
products and evolution equations become ordinary first-order differential equations. By con- 
trast, we are interested in the evolution of truncated moments, defined for a generic function 
f{x) by 

/„(a;o)= C dxx^-^f{x). (2.3) 

JXQ 

The corresponding evolution equations in the nonsinglet case were derived in ref. |]T2|, which 
we will follow in presenting the generalization to the singlet case. One finds immediately that 
the truncated moments of g(x, Q"^) obey the equation 

|g„(xo, Q') = £ dyy-'G^^ qiy, Q') , (2.4) 

where ^ 

Gn{x) = f dzz'^-^Piz) (2.5) 

Jx 
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is perturbatively calculable as a power series in a^. 
Expanding Gn{xo/y) in powers of y around y = I, 



Gr. 



one obtains 



at 



—G f^' 
(9?/P " V 1/ / 



27r 



(2.6) 



(2.7) 



p=0 fc=0 

The key step in the derivation of eq. (p.7| ) is the term-by-term integration of the series expansion. 
This is allowed, despite the fact that the radius of convergence of the series in eq. (^1^ ) is 1 — Xq, 
because the singularity of Gn{xo/y) at y = Xq is integrable (this can be proven |]T3| using the 
Lebesgue definition of the integral). One can then express each power of {y — 1) using the 
binomial expansion, which leads to eq. ( p. 71) . 

Equation (|2.7| ) expresses the fact that, while full moments of parton distributions evolve 
independently of each other, truncated moments obey a system of coupled evolution equations. 
In particular, the evolution of the n*'' moment is determined by all the moments qj, with j > n. 
In practice, the expansion in eq. (p.6|) , because of its convergence, can be truncated to a finite 
order p = M. The error associated with this procedure will be discussed in sect. ^. In this 
case, eq. (p.7|) can be rewritten as 



^qnixo,Q^) 



271 



M 



Y.^nk\xo)qn+k{Xo,Q^) 



{2.i 



k=0 



where 



M 
p=k 



-l)^+%"(xo) 

k\{p-k)\ 



(2.9) 



To solve the system of equations (|2.8| ), it is necessary to include a decreasing number of terms 
(M, M — 1, and so on) in the evolution equations for higher moments (n + 1, n + 2, ...), 
obtaining M + 1 equations for the M + 1 truncated moments . . . , qn+Ai}- We will see in 
the next section that this approximation is fully justified. In this case, the coupled system of 
evolution equations takes the form 



d '^^^ 

-rQk= y2 CkiQi; n<k<n + M, 



(2.10) 



with 



T = dt a[t ) ; a[t) =- 

Jto 



2% 



(2.11) 



where C is now a triangular matrix: 




(M-k+n) 
(^kJ-k 



{I > k) 
\l < k) 



(2.12) 
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In the nonsinglet case, discussed in ref. [0, the matrix elements Cki are just numbers, and 
the matrix C in eq. ( |2.12| ) is triangular, which makes it easy to solve eq. ( p.lO| ) perturbatively. 
In the singlet case, all the steps leading to eq. (|2.12|) are formally the same, but now each 
entry Cm is a 2 x 2 matrix. As a consequence, the matrix C, which is given in terms of partial 
moments of the evolution kernels, is no longer triangular, but has nonvanishing 2x2 blocks 
along the diagonal. This problem can be circumvented, by writing the perturbative expansion 
of C as 

C = Co + aCi + ... ={Ao + Bo) + a{A, + 5i) + . . . , (2.13) 



where A = Aq + aAi is block-diagonal, with 2x2 blocks on its diagonal, 

^kl = Ckk^kl , 



(2.14) 



while B = Bq + aBi, considered as a matrix of 2 x 2 blocks, is upper-triangular with vanishing 
diagonal entries. Now one can define a matrix S that diagonalizes Aq, 



SAnS-^ 



diag(7i, . . . ,72m) • 



(2.15) 



Clearly, S is r-independent, block-diagonal, and easily computed. Equation (|2.10|) can then be 
rewritten as 

^q = Tq, (2.16) 
ar 

where q = S q and T = SCS^^. 

The new evolution matrix T is triangular at leading order (with the same eigenvalues as Aq). 
This is enough to solve the evolution equation to next-to-leading order, as in the nonsinglet 
case of ref. |jl2|. The general solution is worked out in detail in Appendix A; the result is 



g(r) = f/(T,r)g(0) 



(2.17) 



where 



Hk 



->kl 



'ajoy 



ll/bo 



(2.18) 



+ 



kl 



lk-li + bo 



a(0) 



'aioy 



Ik /bo 



- air 



'aioy 



li/bo 



Ij ■ 



In eqs. ( p.l7| , ^TT8D , T is expanded as T = Tq + aTi; R is the matrix which diagonalizes Tq, 
RTqR-^ = diag(7i, . . .72m); finally, Ti = RTiR'^. 

The matrix R can be computed recursively, using the technique applied in ref. and 
proven in Appendix B. One finds 



R 



Rij" 



— X! -^0 



1 ^ 

]] It p=i+l 



(2.19) 
(2.20) 
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which, together with the conditions Ra = 1 and Rij = when i > j, determine the matrix R 
completely. 

The general solution for the parton distributions is then 

g(r) = f/(C,T)g(0), (2.21) 

where 

U{C,t) = S~^U{T,t)S . (2.22) 

The splittting functions and partial moment integrals which should be used in eq. (|2.18|) in 
order compute this solution explicitly are listed in Appendices C and D. 

For the sake of completeness, we describe a different method to solve eq. (|2.8| ). It is imme- 
diate to check that the matrix 

U{C, t)=I+Y. dn... dr^Cin) . . . C(r„) (2.23) 
obeys the differential equation 

±U{C,t) = CU{C,t) , (2.24) 

with the initial condition U{C, 0) = /. In general, eq. ( p.23| ) is not very useful, since it involves 
an infinite sum. In the present case, however, the infinite sum collapses to a finite sum. To 
see this, consider again the decomposition C = A + B, where A is block-diagonal and B is 
upper-triangular. It is easy to prove that 

U{C,T) = U{A,r)U{B,T), (2.25) 

where 

B = U-\A, t)BU{A, t) . (2.26) 

Since A is block diagonal, U{A, r) is also block-diagonal, and it can be computed perturbatively 
using the procedure described in Appendix A. Furthermore, once U{A) is known, the upper- 
triangular matrix B can be computed through eq. ( |2.26|) . Now one can use the fact that upper- 
triangular matrices have the property that their M-th power vanishes. Hence, the solution can 
be expressed as the finite sum 

UiB,T) = I+Y. / dn... dTnB{n)...B{Tn), (2.27) 



n=l 







and from the knowledge of U{B) and U{A) one can determine the solution to the evolution 
equations explicitly. 
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3 Numerical methods and their accuracy 



In this section we will assess the accuracy of our method when the series of contributions to the 



right-hand side of the evolution equation (|2.7| ) is approximated by retaining a finite number 
M of terms. The loss of accuracy due to this truncation is the price to pay for eliminating the 
dependence on parton parametrizations and extrapolations in the unmeasured region. However, 
unlike the latter uncertainties, which are difficult to estimate, the truncation uncertainty can 
be simply assessed by studying the convergence of the series. A reasonable goal, suitable for 
state-of-the-art phenomenology, is to reproduce the evolution equations to about 5% accuracy: 
indeed, we expect the uncertainties related to the parametrization of parton distributions in 
the conventional approach to be somewhat larger (~ 10%)[|. Notice that there is no obstacle to 
achieve a higher level of precision when necessary, by simply including more terms in the relevant 
expansions. To this level of accuracy it is enough to study the behavior of the leading order 
contribution to the evolution equation: indeed, next-to-leading corrections to the anomalous 
dimension are themselves of order 10%. We have verified explicitly that the inclusion of the 
next-to-leading corrections does not affect our conclusions. 

We can compare the exact evolution equation ( |2.7| ) with its approximate form, eq. ( p.8[ ), by 
defining the percentage error on the right-hand side of the evolution equations for the quark 
nonsinglet, singlet and gluon: 



■TJNS 



■n,M 




(3.1) 



(3.2) 



(3.3) 



where N'Ns,T.,g are the exact right-hand sides of the evolution equation (PT7|). We study the 
dependence of the percentage error on the value of M for typical values of the cutoff Xq and 
for representative choices of test parton distributions. In particular, we parametrize parton 
distributions as 



Nx~ 



[l-x 



,/3 



(3.4) 



We begin by choosing, as a representative case, /3 = 4 and a = 1 for the singlet distributions 
and a = for the nonsinglet. The nonsinglet is assumed to behave qualitatively as qNs ~ 
xg ~ xS, in accordance with the behavior of the respective splitting functions. Furthermore, 



^Notice that this is not the uncertainty associated with evolution of a given parametrization with, say, an 
x-space code; rather, it is the uncertainty associated with the choice of the parametrization, and with the bias 
it introduces in the shape of the distribution. 
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the normalization factors for the singlet and gluon are fixed by requiring that the second 
moments of Tj{x,Q^) and g{x,Q'^) are in the ratio 0.6/0.4, which is the approximate relative 
size of the quark and gluon momentum fractions at a scale of a few GeV^. We will then show 
that changing the values of a and (3 within a physically reasonable range does not affect the 
qualitative features of our results. 

The accuracy of the truncation of the evolution equation is determined by the convergence of 
the expansion in eq. ( |2.6| ). Because this expansion is centered at ?/ = 1, and diverges at y = Xq, 
the small y region of the integration range in eq. ( p.4|) is poorly reproduced by the expansion. 
Hence, even though the series in eq. (|2.7|) converges, as discussed in sect. 2, the convergence will 
be slower for low moments, which receive a larger contribution from the region of integration 
y ^ xq. In fact, for low enough values of n, the convolution integral on the right-hand side 
of the evolution equation (|2.4|) does not exist: this happens for the same value for which the 
full moment of the structure function does not exist, i.e. n < 1 in the unpolarized singlet and 
< in the unpolarized nonsinglet and in the polarized case. Therefore, we concentrate on 
the lowest existing integer moments of unpolarized distributions, i.e. the cases n = 2, 3 for the 
singlet distributions, and correspondingly n = 1, 2 for the nonsinglet, which are the cases in 
which the accuracy of the truncation will be worse. 

The values of 7^^"^^'^, computed at leading order with xq = 0.1, are shown in table |l]. The 
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0.63 


0.43 


0.55 


0.16 


0.12 


0.16 
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0.49 


0.36 


0.38 


0.13 


0.10 


0.12 


20 


0.34 


0.27 


0.26 


0.10 


0.08 


0.08 


40 


0.20 


0.17 


0.17 


0.06 


0.05 


0.05 


70 


0.12 


0.10 


0.10 


0.04 


0.03 


0.03 


100 


0.09 


0.07 


0.07 


0.03 


0.02 


0.02 


150 


0.06 


0.05 


0.05 


0.02 


0.01 


0.01 



Table 1: Values ofTiJ^^^'^ for xq = 0.1 and different values of n and M . 

table shows that nonsinglet moments of order n behave as singlet moments of order n — 1. 
This is a consequence of the fact that, as discussed above, the convergence of the expansion 
is determined by the singularity of the integrand (?„ {xo/y)q{y) of eq. ( P^ ) as y ^ Xq; near 
y = Xq, the function (?„ (xo/y) is well approximated by the singular contribution log (1 — xo/y), 
while parton distributions carry an extra power of y~^ in the singlet case in comparison to the 
nonsinglet. We also observe in table 1 that, as expected, the convergence is slower for the lowest 
moments, and rapidly improves as the order of the moment increases. This rapid improvement 
is a consequence of the fact that the convergence of the expansion of G{xQ/y) is only slow in 
the immediate vicinity of the point y = xq, and the contribution of this region to the n-th 
moment is suppressed by a factor of Xg"^. Due to this fast improvement, the approximation 
introduced by including one less term in the expansion as the order of the moment is increased 
by one, which is necessary to obtain the closed system of evolution equations ( p.lO| ), is certainly 
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justified. 

The 5% accuracy goal which we set to ourselves requires the inclusion of more than 100 
terms for the lowest moment, but only about 40 terms for the next-to- lowest. The computation 
of series with such a large numer of contributions does not present any problem, since the 
splitting functions are known and their truncated moments are easily determined numerically. 
The implications of this requirement for phenomenology will be discussed in the next section. 

We can now study the dependence of these results on the value of the truncation point Xq 
by plotting the exact and approximate right-hand side of the evolution equations as a function 
of Xq, as shown in fig. |I]. The figures show that the case Xq = 0.1 studied in table 1 is a generic 
one between the two limiting (and physically uninteresting) cases Xq = and Xq = 1, where the 
approximation is exact. In fact, with this particular choice of parton distributions, xq = 0.1 is 
essentially a worst case and the error estimates of table 1 are therefore conservative. 

An interesting feature of these plots is the presence of zeroes of the lowest moment evolution 
at Xo = in the nonsinglet and around xq ~ 10~^ in the gluon case. The physical origin of these 
zeroes is clear. At leading order, the first nonsinglet full moment does not evolve. On the other 
hand, the second gluon full moment grows with Q^, while higher gluon full moments decrease, 
i.e. the gluon distribution decreases at large x; this implies that the second truncated moment 
of the gluon must decrease for a high enough value of the cutoff Xq, while it must increase for 
very small Xq; its derivative is thus bound to vanish at some intermediate point. Of course, 
the phenomenology of scaling violations (such as a determination of as) cannot be performed 
at or close to these zeroes, where there is no evolution. From the point of view of a truncated 
moment analysis, this means that the value of Xq should be chosen with care in order to avoid 
these regions. 

Finally, in table |^ we study the dependence of our results on the form of the parton dis- 
tributions, by varying the parameters a and (3 within a reasonable range. Of course, parton 
distributions which are more concentrated at small y give rise to slower convergence. However, 
we can safely conclude that the effect of varying the shape of parton distributions is generally 
rather small. We have also verified that varying the relative normalization of the quark and 
gluon distributions has a negligible effect on the convergence of the series, even though it may 
change by a moderate amount the position of the zeroes in gluon evolution discussed above. 



4 Techniques for phenomenological applications 

So far we have discussed scaling violations of parton distributions. In a generic factorization 
scheme, the measured structure functions are convolutions of parton distributions and coeffi- 
cient functions. When taking moments, convolutions turn into ordinary products and moments 
of coefficient functions are identified with Wilson coefficients. In the present case, however, as 
shown in eqs. (|2.4| - |2.7|) , truncated moments turn convolutions into products of triangular ma- 



trices. Hence, in a generic factorization scheme, truncated moments of parton distributions are 
related to truncated moments of structure functions by a further triangular matrix of truncated 
moments of coefficient functions. 

This complication can be avoided by working in a parton scheme 0], where the quark 
distribution is identified with the structure function F2. This still does not fix the factorization 
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Figure 1: Right-hand sides of the evolution equations for the first and second truncated moments 
of the nonsinglet distribution, and for the second and third moments of singlet distributions. 
The overall scale is set by as{2GeV^). 
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0.01 


1.5 


4.0 


0.12 


0.05 


0.11 
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0.05 


0.02 


0.05 


0.02 


1.5 


6.0 


0.16 


0.07 


0.15 


0.06 


1.0 


6.0 


0.12 


0.05 


0.12 


0.05 


0.5 


6.0 


0.09 


0.03 


0.08 


0.03 



Table 2: Values o/7^^^^ and 'RPn m f^"^ = 0.1 and different choices of parton distribution 
parameters. 



scheme completely in the gluon sector. One way of fixing it is to use a "physical" scheme, where 
all parton distributions are identified with physical observables |^. This may eventually prove 
the most convenient choice for the sake of precision phenomenology, once accurate data on 
all the relevant physical observables are available. At present, however, the gluon distribution 
is mostly determined from scaling violations of F2, so within the parton family of schemes 
the choice of gluon factorization is immaterial. Here we will fix the scheme by assuming that 
all moments satisfy the relations between the parton-scheme gluon and the MS quark and 



gluons imposed by momentum conservation on second moments [|I4|. This is the prescription 
used in common parton sets, and usually referred to as DIS scheme. Explicit expressions of 
Altarelli-Parisi kernels in the DIS scheme are given in Appendix |C[ 

With this prescription, the phenomenology of scaling violations can be studied by computing 
a sufficiently large number of truncated moments of structure functions, so as to guarantee the 
required accuracy. If the aim is, for instance, a determination of from nonsinglet scaling vio- 
lations, all we need is a large enough number of truncated moments of the nonsinglet structure 
function. Once an interpolation of the data in the measured region is available, the determi- 
nation of such truncated moments is straightforward. This interpolation can be performed in 
an unbiased way using neural networks, as already mentioned in the introduction. One may 
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wonder, however, whether the need to use the values of very high moments wouldn't be a prob- 
lem. Indeed, very high moments depend strongly on the behavior of the structure function at 
large y, which is experimentally known very poorly. Furthermore, it seems contradictory that 
scaling violations of the lowest moments should be most dependent on the structure function 
at large y. 

This dependence is only apparent, however. Indeed, eq. (|2.7|) for, say, the nonsinglet first 
moment can be rewritten as 



The need to include high orders in the expansion in eq. (|2.^ ) is due to the slow convergence of 
the series in eq. ([4.1|) , in turn determined by the fact that Gn{xo/y) diverges logarithmically at 
y = Xq. Correspondingly, the right-hand side of the evolution equation depends significantly on 

for large values of p, which signal a sensitivity to the value of q{y, Q^) in the neighborhood 
of the point y = xq. The dependence on high truncated moments g„ is introduced when 
is re-expressed in terms of g„, by expanding the binomial series for {y — ly. Since this re- 
expansion is exact, it cannot introduce a dependence on the large y region which is not there 
in the original expression. The high orders of the expansion do instead introduce a significant 
dependence on the value of the structure function in the neighborhood of xq, which can be 
kept under control provided xq is not too small, i.e. well into the measured region. There 
is therefore no obstacle even in practice in performing an accurate determination of as from 
scaling violations of truncated moments. 

Let us now consider a second typical application of our method, namely the determination 
of truncated moments of the gluon distribution. In particular, the physically interesting case 
is the lowest integer moment, i.e. the momentum fraction in the unpolarized case or the spin 
fraction in the polarized case. The need to include a large number of terms in the expansion 
of the evolution equations seems to imply the need to introduce an equally large number of 
parameters, one for each gluon truncated moment. This would be problematic since it is appears 
unrealistic to fit a very large number of parameters of the gluon from currently available data 
on scaling violations. We may, however, take advantage of the fact that the dependence on high 
order truncated moments is fictitious, as we have just seen, and it rather indicates an enhanced 
sensitivity to the value of q{y) as y ^ Xq. This suggests that a natural set of parameters to 
describe the gluon distribution should include the first several truncated moments, as well as 
further information on the behavior of the distribution around the truncation point xq, such as 
the value of the distribution (and possibly of some of its derivatives) at the point xq- 

To understand how such a parametrization might work, notice that if q{y) is regular around 
y = Xq, then it is easy to prove that 




(4.1) 



where 




(4.2) 



p—fOO 



lim 



/4 dyjy - lyqjy) 
qi^o) l'dy{y-l)P 



1 , 



(4.3) 
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by Taylor expanding q{y) about y = 
appears on the r.h.s. of eq. ( |4.1| ) by 



We may therefore approximate the series which 



S{xo,no) 



^ /' dy{y - ir<iiy) + E 



p=0 



p=no 



q{xo) / dy{y 



(4.4) 



Equation (|4.4| ) describes the evolution of the first truncated moment of q{y) in terms of the 
first no truncated moments and of the value of q{y) at the truncation point xq. Of course, the 
approximation gets better if uq increases. It is easy to check that when xq = 0.1 the accuracy 
is already better than 10% when uq ~ 7. This means that a parametrization of the distribution 
in terms of less than ten parameters is fully adequate. It is easy to convince oneself that 
this estimate is reliable, and essentially independent of the shape of the distribution q{y). In 
fact, because slow convergence arises due to the logarithmic singularity in Gn{xo/y), we can 
estimate the error of the approximation in eq. ([4.4|) by replacing the functions gp{xo)/p\ with 
the coefficients of the Taylor expansion of log(l — Xo/y) in powers ofy — 1, which we may denote 
by gl{xo)/p^-- The error is then 



p=no 



p\ 



dy{y-lf {q{y) - gixo)) 



XQ 



(4.5) 



< I'dy 

IXQ 



iog(i-Wz/)-"E 



\q{y) - g(xo)| . 



The expression inside the first absolute value on the r.h.s. of eq. ( [4. 51) is just the error made 
in approximating the logarithm with its Taylor expansion around y = 1\ thus, it is a slowly 
decreasing function of no, it is integrable, and the integral receives the largest contribution 
from the region ?/ ~ xo; the second absolute value, on the other hand, is a bounded function 
of y in the range Xq < y < 1, which vanishes as ?/ Xq for any choice of q{y). These two 
facts combine to limit the size of the error. One can check directly that, choosing for example 
qiy) = (1 — yY as in the previous section, the accuracy is better than 10% with no ~ 10 and 
Xo = 0.1, in agreement with the previous estimate. One may also verify that, as expected, 
changing the shape of qiy) does not significantly affect the result. 

We conclude that there is no difficulty in using the evolution of truncated moments for a 
direct extraction of the lowest truncated moment of the gluon distribution, provided only the 
higher moments of the distributions, which are auxiliary quantities needed in the extraction, 
are parametrized in an effective way. We have seen that this is possible with a reasonably small 
number of parameters. 



5 Outlook 

In this paper, we have discussed the solution of the Altarelli-Parisi evolution equations of per- 
turbative QCD by projecting parton distributions on a basis of truncated moments. We have 
seen that truncated moments give us a compromise between standard Mellin moments, which 
satisfy simple linear evolution equations, but are not measurable, and parton distributions 
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themselves, which are measurable, but satisfy integro-differential evolution equations. In this 
respect, projecting on a basis of truncated moments is akin to projecting on a basis of orthogonal 
polynomials but with the added advantage that truncated moments are physical observ- 
ables. We have further shown that evolution equations for truncated moments can be solved 
to arbitrarily high accuracy by using a sufficiently large basis of truncated moments, and we 
have discussed how this formalism can be exploited to perform a model-independent, unbiased 
analysis of scaling violations, which does not rely on parton parametrizations. We have also 
collected all technical tools which are needed for this analysis, and discussed the reliability of 
the approximations which are necessary in order to implement it in practice. Such an analysis 
could be used for the determination of as and for the direct measurement of the contribution to 
the moments of the gluon distribution from the experimentally accessible region. In both cases, 
it would offer the advantage of allowing reliable estimates of the uncertainty on the result. In 
order to be effective, these phenomenological applications will need an unbiased interpolation 
of the data, such as could be achieved by means of neural networks. Phenomenological stud- 
ies along these lines are currently under investigation and will be presented in forthcoming 
publications. 
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A Solution of matrix evolution equations at NLO 

In this Appendix, we find the general solution of the equation 

A, = C,, (A.1) 

where g is a vector with M components, and C is a generic M x M matrix. The usual QCD 
evolution equations are special cases of this equation in which M < 2. We will assume that C 
has a perturbative expansion in powers of a parameter a(r): 

C = Co + a(r)Ci + ... , (A.2) 

with 

^ -6oa (l + 6ia +...) . (A.3) 



dT 

For QCD applications 



as 1 



a = ^; ^ = ^ / dt'asit') , (A.4) 



to 



with t = log((5^/A^^^), and 

6o_33-2n/_ bi_ 153 - 19n/ 



2tt Utt ' 2n 27r(33 - 2n 



f) 



(A.5) 
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The solution of eq. ( |A.1|) can be obtained perturbatively. Expanding q to order a, 



we find 



q = qo + aqi 



— go = C'o Qo 
dr 



d 

dr 



qi = (Co + bo) qi + Cigo 



(A.6) 

(A.7) 
(A.8) 

(A.9) 

gi(r) = i?-ie(^+''o)"i?gi(0) + i^-ig^^+^o)^ T rfa e'^^+^'o^'^C'ie^'^ i?go(0) , (A. 10) 

Jo 

where the matrix R diagonahzes Co, 

RCoR-^ = diag(7i, . . . , 7m) = 7 , (A.ll) 

and 



The solutions of eqs. (|A.71 , |A.8| ) are 

go(r) = /?-V^/2go(0) 



Ci = RCiR ^ . 



(A.12) 



Collecting these results, and noting that aexp(6or) = a(0), up to terms of order a^, we can 
write the solution as 



g(r) = U(C,T)q(0) = R'^ e^" + ae(^+^°)" / c/a e-(^+^")"C'ie^" 

L Jo 

with the initial condition 

g(0) =go(0) + a(0)gi(0). 
The explicit expression of U{C,t) is 



Rq{0) 



which, expanded to next-to-leading order reduces to 



o-r^T _ „(7™+bo)r 



7" — 7™ — Bq 



R. 



(A.13) 
(A.14) 

(A.15) 



Uij{C,T) = ^ 



+ 



a[T)j 
7"^ - 7" 6o 



a(0) 



'a(O)' 



, a r 



Im/bfj 



^a(O)^ 



In /bo 



>R 



(A.16) 



In the case of standard QCD evolution equations the matrix Cq is at most 2x2 and is eas- 
ily diagonalized. In the cases treated in this paper, the matrix Co is triangular and can be 
diagonalized using the methods discussed in the next Appendix. 
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B Diagonalization of triangular matrices 



In this Appendix, we show how to construct the matrix R which diagonahzes a generic n x n 



triangular matrix T by means of the recursion relations eqs. (|2.19| , |2l20| ). The matrix R is 
defined by the requirement that 

i?Ti?-i = diag(7i,...,7n) , (B.l) 

where the matrix T is upper triangular, i.e. Tij = if i > j. It is easy to see, by solving the 
secular equation, that the eigenvalues 7^ of T coincide with its diagonal elements, 

7i = Tu . (B.2) 

Now, define eigenvectors associated to the j-th eigenvalue Tjj, with components Vi^: 

n 

J2 TikVk^ = jjVi^ . (B.3) 

k=l 

Clearly, the matrix R^^ coincides with the matrix of right eigenvectors, {R~^)ij = Vi\ while 
the matrix R coincides with the matrix of left eigenvectors J2k=iV-' kTki = '^jV^'i, Rij = v^j- 
The eigenvector condition eq. ( [B.3| ) immediately implies that the j-th component of the j-th 
eigenvector is equal to one: Vj^ = 1. Furthermore, it is clear that eq. (|B.3|) can only be satisfied 
if all components Vk-' of the j-th eigenvector with k > j vanish, 

v/ = 1 ; Vk^ = if k>j. (B.4) 

Using eq. ( p.4|) and the fact that the matrix T is triangular, eq. (p.3| ) can be written as 

j 

J^TikVk^ = 'jjVi^ . (B.5) 

k=i 

Substituting the explicit form of the eigenvalues, eq. ( [B.2| ), and identifying Vi^ = {R~^)ij, this 
is immediately seen to coincide with eq. ( 2.201) . Furthermore, using the condition vj^ = 1, this 



equation can be viewed as a recursion relation which allows the determination of the {k — l)-th 
element of once the k-th element is known, which is what we set out to prove. The same 
argument, applied to the left eigenvectors, leads to the expression in eq. ( p.l9|) for R. 



C Splitting Functions in the DIS scheme 

In this Appendix, we give the explicit expressions of the Altarelli-Parisi kernels in the DIS 
scheme We define the non-singlet splitting functions as 

P±(x) = P,,{x) ± P,g{x) = Pf (x) + aPi'\x) + ... (C.l) 
and the singlet 2x2 matrix of splitting functions as 

P{x) = P(°) (x) + aP(i) (x) + . . . . (C.2) 
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The MS LO and NLO kernels are given for example in eqs. (4.94) and (4.102)-(4.112) of ref. [||, 
whose notation and conventions we have followed throughout this paper. The splitting functions 
in the DIS scheme can be constructed from these by a change of factorization scheme. 

To next-to-leading order, a generic change of factorization scheme for the splitting functions 

is 



(C.3) 
(C.4) 



where E and E^^ are functions of a;, and the commutator in eq. ( |C.4| ) is defined in terms of 
convolutions, as 



P(°) = E® P(°' - PC^) ® E 



(C.5) 



and is thus in general nontrivial to compute. 

The transformation that takes from the MS to the DIS scheme, defined as in ref. |n 
given by 



IS 



E' 



E{x) 



(C.6) 
(C.7) 



C"?(x) 2nfC3{x) 
-Ci^x) -2nfC3{x) 

where C"^(x) and C^{x) are the next-to-leading terms of the quark and gluon coefficient functions 
for the unpolarized deep- inelastic structure function E2{x,Q^) in the MS scheme: 



x{e) — 
Jx y 



X 



where (e^) = J2i=q,qej /{2nf). Exphcitly, 

/ logf 1 — x) 



X 



5 1 - - + aCM - 



j:{y,Q') + 2nfaC^i-\g{y,Q' 



/•I dy 








Jx y 




^ y) 





QNsiy,Q^) + ■■■ 



Cix) 



1 



X 



X 



X 



X 



log X + 3 + 2x 



IT 

y ^ 2 



— (1 + x) log(l — x) 
9' 



(5(1 -x) 



C^x) 



T, 



R 



f (1 - x^) + x"^) log -Sx"^ + 8x-l 

^ ' X 



(C.9) 
(C.IO) 



Notice that the expression for the scheme change given in eq. (3.9) of ref. lacks a factor of 
2nj in the qg entry because it refers to a single flavor. 
The explicit expression of the commutator is 

r2\ 



E,PW 



11 



Cpnf 



1 / vr' 
-3Li2(a;) (l + 4x-4x^) log^x 

2 \ 3 



'4 5 ^ 

h7x 

3x 2 



22x^\ , , X / 3 22x2 

log(l -x)- - + 9x 

3 / ' \2 3 



logx 



;i-2x + 2x2) log(l-x)log- 



-X 2 55 20 X 
x^ 3x 6 3 



IX' 



,(C.ll) 
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19 



/ 1 TT^ 

(l + 4x + 4x2) (Li2(x) + -log'x- — 



{2 + Ax-6x'^) log(l-a;) + (2 + 16a;+10a;^) logo; 



11 , 59a;2 



+ 

+ CpUf 



-(l-2x + 2x2) 2Li2(a;) +log^(l -x) - 



7r 



(1 + lOx- lOx^) log{l-x)-{l-4:X+lQx^) loga;-6 + llx-8x2 



+ 2C^n/ 



- (1 + 4a;) (Li2{x) + ^ log^ a; 



2 1 79a;^\ , 



+ (1 -2a; + 2x^) log(l -x)log 



1 — a; 



X 



'1 ... 3ixn , 

- + 16a; — I loga; 



1 43 121a; 



^ 3a; 12 



+ 7r^x + 



281x^ _ Tr^a;^ 
12 3~ 



(C.12) 



boC'^{x) + CFnf 



(l + x) (log^a; + 2Li2(a;)) 



+ ^(^20 - - ^ - 2x - n^x - 16x^j 

/ . 4a;2\ , / 4 4a;2 



log(l - a;) 

3 1 2 
+ - log X 



6 5a; 



12-- + :^llog(l-.)+(2-2-,],og^(l-,)-^_^ + | 



- 2 2 X 

\ X 



- (l + a;2) 



^ loga;log(l — a;) — ^1 2a;^ ^i2{x) + + "^-^^ loga; 

log^(l — a;) ^ I log^ X ^ log a; log(l — a;) 



1 -a; 



1 -a; 



1 -a; 



9, loga; 



27 2 7r2 27rV 



1 — X 



+ 



'^log(l — x) 

1 — X 
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+ 2|^l---2xjLi2(a;)- 

+ 2x 9x + ^— 

1 — X \ 3 



logxlog(l -x) 

iX : h OX 



1 — X 



log (1 — x) 

1 — X 



1 — X 



;i + 2x^)(3 logx + log^x) 



+ [- - 2xA log^(l - s) + (2 - - - 2x + 4a;^ j loga;log(l - x) 



X 



6 



+ 17 4a; + 6a;^ log(l - x) - 2 - 6x + — — + 8a; + 



2n^x^ 



X 



(C.13) 



Finally 



Note that 



99 



(C.14) 



99 



denotes the qg matrix element of the commutator; i.e., using the conven- 



tions of ref. H], upon scheme change 



2nyP« ^ 2n;P« + k P(°) 



99 



boEqg . 



(C.15) 



D Truncated moment integrals 

The integrals which are needed in order to compute the Mellin moments of the NLO Altarelli- 
Parisi splitting functions are well known [|l^]. Here we list all the truncated moment integrals 
which are necessary in order to determine the evolution kernels Gn{x), eq. ( |2.5| ), from the 
expressions of the splitting functions given in appendix C. The triangular anomalous dimension 
matrices C^i, eq. ( p^.l2D , can be easily determined from these formulas by Taylor expansion. 



dz z 



n-l 



dz z^ ^ log z 
r dzz'^-Hog^z 

Jx 

dz z'^-Hodl - z) 



dzz^ Mog^logfl — 2) 



1 - z 



1 — X^' 



n 



1 — X" (1 — 72 logx) 



2 — x" (2 — 2 n log X + log^ x) 



n+l 



(D.l) 
(D.2) 
(D.3) 



;Pi(n + l,l;n + 2;x) 



n (n + 1) 

+ (n + 1) (-iE + log(l - x) + i)^^\n + 1; 



(D.4) 



-iE-n x"+^ <l>(x, 2, n + 1) - x" log(l - x) (D.5) 



-l-rax" log(l — x) logx + x""*"^ $(x, 1, n + 1) {n logx — 1) 
+^(o)(^) + n - n [V^(°)(n + 1)]^ + nV^«(n + 1) 



X 



+ x"+^$(x,2,n + 11 



X 



n 



logx 
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1 — z 



/ dzz 

Jx 



1-Z 



fdzz'^-^-^ 

Jx 1 + Z 



L 



l + z 



Cdzz-^'^ 

Jx 1 + Z 



"dzz--^^''^-'^ 



l + z 



I dzz^ ^ log z\og{l + z) 

Jx 



-X 



1, n + 1) log X - -0^^^ (n) , 



_2x"+i$(x,3,n + l) + ^ logx 



+2 x''^^ <l>(x, 2, n + 1) log X log^ X 

n 

1, n + 1) log^ X - -0^^^ (n) , 

r'dzz^-'(-^) = -7E-— 2Fi(l,n;n+l;a;)-^W(n) , 
\1 — 2;/+ n 



Li3(l - x) - log(l - x)Li2(l - x) 

n-l „i 

ci2;2;*^~^ log2;log(l - , 



fe=l 



1 



2Fi(n, l;n+ 1; -x) 



2 



n+ 1 



C(3) logxlog^(l + a;) log^(l + x) 



8 



+ 



+ log(l + x)Li2 (y^) + Li3(-2;) + Lis ( 

n-l „i 

-^(-i)^-W (i^z^-Mog^iog(i + ^) 



1 + X 



fe=l 

2 a;" 



+ 2a;"+i$(-a;,3,n + l) 



2 

9 



x" log X - 2 x"+^ ^(-a:;, 2, n + 1) log x 



log X + <l>(-x, 1, n + 1) log X 

n 



(2) 



(2) 



n + 1 



= (-1) 



n-l 



I^-^log2 + log(l + x)Li2(-x) 



+ logx log2(l + 2;) + y log(l + x) + 2Li3 (y^) 



dzz^ ^Liof-^) 



fc=i 



4 n 0;"+^ $(-a;, 2, n + 1) - 4 log 2 



(D.6) 



(D.7) 
(D.8) 



(D.9) 



(D.IO) 



(D.ll) 



(D.12) 



(D.13) 



+4x''+^$(-x, + 1) (-1 + n logx) +4x" log(l + x) 
-4nx" logx log(l + x) + 2^(°) + -^V"^"^ 



n + 1 
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1-z 



+ 



Jx 



dzz^'-^Uoi-z) 



f 

Jx 



dz Z''-^U2{Z) 

Jx 



1-Z 



+ 
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(D.14) 

(n - 1) (1 - x) [4^3(1, 1, 1, 2 - n; 2, 2, 2; 1 - x) (D.15) 

-3^2(1, 1, 2 - n; 2, 2; 1 - x) log(l - x)] + , 

^ x''+\Fi{n + l,l;n + 2;-x) 



n2 (n + 1) 
-in +1)1 
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log 2 + a;" log(l + 



(D.16) 
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n+ 1 



nx"' Li2(— x) 



z"-' log^(l -z) = 2 (1 - x) 4^^3(1, 1,1,1- n; 2, 2, 2; 1 - x) 



-3F2(l,l,l-n;2,2;,l-a;)log(l-x) 
+^^\og\l-x), 



n 
1 
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log(l - xY 



n-2 



(i-^)E 

3=0 L 



2 4F3(l,l,l,-j;2,2,2;l-a;) 



■2 3F2(l,l,-j;2,2;l-x) log(l - x) 

l-x^+^ . 2. 



(j + l)(l-x) 



log (1 — x) 



(D.19) 



00 



A;=l 



= '''"'lil"';.'"' ; -TE = « 0.677216 ; 



^^^(ai, . . . , Op; 61, ■ • ■ , 6g; 2;) = E 



(/c + ai) . . . (/c + ttp) z'^ 



^0 (/c + 61) • • • (/c + 6g) k\ ■ 



(D.20) 

(D.21) 
(D.22) 

(D.23) 
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Note that the integrals in eqs. (p.9| ), ( p. 11 ), ( p.l3|) and (p.l9|) are vahd only when n is a 
positive integer, which is what is usually needed. If an expression for real or complex n is 
necessary, one should compute these integrals numerically. All other integrals are given in a 
form that immediately generalizes to complex n. The special functions in eqs. ( p.20[ - [IX23D are 
available in algebraic manipulation programs of common use. 
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